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The effects of passive scalar anisotropy on subgrid-scale (SGS) physics and modeling 
for Large-Eddy Simulations are studied experimentally. Measurements are performed 
across a moderate Reynolds number wake flow generated by a heated cylinder, using 
an array of four X-wire and four cold-wire probes. By varying the separation distance 
among probes in the array, we obtain filtered and subgrid quantities at three different 
filter sizes. We compute several terms that comprise the subgrid dissipation tensor of 
kinetic energy and scalar-variance and test for isotropic behavior, as a function of filter 
scale. We find that whereas the kinetic energy dissipation tensor tends towards isotropy 
at small scales, the SGS scalar-variance dissipation remains anisotropic independent of 
filter scale. The eddy-diffusion model predicts isotropic behavior, whereas the nonlinear 
(or tensor eddy diffusivity) model reproduces the correct trends, but overestimates the 
level of scalar dissipation anisotropy. These results provide some support for so-called 
mixed models but raise new questions about the causes of the observed anisotropy. 



1. Introduction 

The statistics and general structure of passive scalars in turbulent flows differ signif- 
icantly from those of the turbulent velocity field. In particular, conclusive experimental 
(e.g. Stewart 1969, Sreenivasan, Antonia & Britz 1979, Mestayer 1982, and Mydlarski 
& Warhaft 1998a) and numerical (e.g. Holzer & Siggia 1994) evidence shows that struc- 
ture functions and the derivative skewness of the scalar field do not follow predictions 
from isotropy at inertial and dissipative scales, in the presence of a mean scalar gradient. 
In particular, the deviations are thought to be related to "ramp and cliff structures" 
and to imply a direct effect of large-scale structures on small-scale structures. The data 
relevant to this question has been reviewed by Sreenivasan (1991) and more recently 
by Warhaft (2000). Moreover, scalar spectra display a distinctly less universal structure 
than velocity spectra. This is manifested both in terms of spectral exponents, as well 
as in terms of the dimensionless spectral coefficient cg. The latter varies from values 
near 0.4 in many experiments in the atmospheric surface layer and grid turbulence (see 
Sreenivasan 1991, Sreenivasan 1996 and Warhaft 2000) to cg ~ 1.8 for other atmospheric 
measurements (Antonia, Ould-Rouis, Anselmet & Zhu 1997, and Antonia, Xu & Zhou 
1999). From a fundamental point of view, these observations challenge the Kolmogorov 
cascade phenomenology for the transfer of scalar-variance from large to small scales. In 
the classical phenomenology, the multiplicity of separate "eddy breakdown" events is 
assumed to gradually uncouple the small from the large scales, allowing the former to 
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tend to a universal and isotropic structure more or less independent of the large scales. 
On the contrary, the scalar-field measurements suggest a direct linkage among largest 
and smallest scales. 

The coupling among scales is an important ingredient in Large-Eddy Simulations 
(LES). In LES, the turbulent fields (velocity and scalar) are decomposed into large and 
small (subgrid-scale, SGS) scale contributions by means of a spatial low-pass filter of 
characteristic width A. The resulting equations which can be numerically discretized 
with mesh-spacing of the order of A require closure of the unresolved momentum fluxes 
(SGS stress tensor, Tij = "UiUj — UiUj, where a tilde represents filtering at scale A) and 
scalar fluxes (e.g. the SGS heat flux, qj = 9uj — 9uj where 6 is the passive scalar field). 
The promise of LES is often predicated upon small-scale universality and isotropy, and 
the absense of a strong coupling across disparate length scales. Hence, the observed 
anisotropy of the scalar field seems to pose a challenge to the very foundation of LES. 
While these deviations from classical phenomenology are now quite well established, little 
is known about their impact on the closure problem for LES. The present work quantifies 
the implications of small-scale scalar anisotropy on quantities that describe subgrid-scale 
physics and directly affect modeling for LES. 

As reviewed in Meneveau & Katz (2000), the most important statistical property of 
the fiuxes r^j and cjj is how they affect the mean kinetic energy and scalar-variance 
budgets of the resolved fields. Specifically, their dominant effect is through the kinetic 
energy and scalar-variance dissipations that arise from interactions between subgrid and 
resolved scales. Therefore, in the present study, we mainly focus on tht^ so-called SGS 
kinetic energy dissipation —{TijSij) (Piomelli et al. , 1991) and scalar- variance dissipation 
— {qjGj) (Porte-Agel, Meneveau & Parlange 1998). Here Sij = \{dui/dxj + duj/dxi) 

and Gj = dO/dxj are the resolved strain-rate tensor and scalar gradient, respectively. 
The SGS dissipation rates represent the flux (cascade) of kinetic energy or scalar- variance 
from resolved towards subgrid scales (when positive). When A pertains to the inertial 
range, and when the flow is in equilibrium, one expects the mean SGS dissipation to 
equal the molecular dissipation rate. 

Deviations from isotropy in the context of SGS dissipation can be probed by measuring 
the isotropy level of the tensors —{TijSmn) and —{qiOj), as a function of scale. The main 
question to be addressed in this work is whether the approach to isotropy (if it exists) 
is the same for kinetic energy and scalar-variance dissipation tensors. Another goal is to 
test the ability of two popular model classes (eddy diffusivity and nonlinear models, see 
Meneveau & Katz 2000 for a review) to reproduce the observations. In order to span a 
sizeable range of inertial-range filter scales, a sufficiently high Reynolds number must be 
considered. Hence, this study is based on experimental data (as opposed to DNS which 
is limited to small Reynolds numbers). The study is performed in a canonical shear flow, 
the heated cylinder wake. 



2. Experiment apparatus and flow characteristics 

Experiments were performed in the retmn type Corrsin Wind Tunnel (Comte-Bellot & 
Corrsin 1966). A heated smooth cylinder of diameter D = 4.83 cm was located horizon- 
tally at the centerline of the test section. The measurement location in the streamwise 
direction (a;i) was fixed at xi/D = 25. To obtain the filtered and SGS quantities, an 
array of four custom-made miniature probes was used. Each probe was composed of one 
X-type hot-wire and one I-type cold-wire for the velocities in the xi — X2 plane and the 
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temperatures, respectively. Here X2 is the "cross-wake" direction, i.e. perpendicular to 
xi and to the cylinder axis. 

The separation distance h between the probes in the cross-wake direction X2 could be 
adjusted manually between 5 and 20 mm. Three configurations (A = 2/i, with /i = 5, 10 
and 20 mm) were used in the present study. A 2.5 fim platinum-coated tungsten wire 
which had been copper-plated was soldered on to the X-wire prong ends and etched, 
yielding an active length-to-diameter ratio of about 200. The wire spacing between the 
hot-wires was 0.5 mm. A 0.625 /im silver-coated pure platinum wire for a cold-wire 
sensor was soldered on the I-wire prong ends and subsequently etched. To minimize the 
low frequency a mplitude attenu ation, the active length-to-diameter ratio was about 1000 
as suggested by Bruun (1995) . The separation distance between the cold-wire and its 
nearest hot-wire was 0.9 mm so that the thermal effect from the hot-wire on the cold- 
wire was negligible. The signals were low-pass filtered at a frequency of 20 kHz and 
sampled at fs = 40 kHz. Sampling time was 60 seconds, so the total number of data 
points per channel for each measurement location was 2.4 x 10^. The array was traversed 
across the wake, and data were recorded at 17 discrete cross-wake locations from the 
centerline to the wake edge at increments of 14.4 mm. 

At the measurement location of Xi/D = 25, the mean centerline velocity [Ucl) was 
13.6 m/s, the defect velocity {Ud — Uoq — Ucl) was 4.4 m/s, the defect temperature {dd = 
GcL — Ooo) was 0.61 °C, and the half-width of the wake was £ — 0.08 m. To get the spatial 
quantities along the streamwise direction from the temporal data, Taylor's hypothesis was 
invoked. The turbulence intensity of the streamwise velocity at the centerline was about 
13.3%. The molecular kinetic energy dissipation at the centerline (ecl) was 87 m^s~^, 
and the molecular scalar-variance dissipation at the centerline {eecL) was 0.65 °C^s~^. 
The latter two variables were obtained from (corrected) third-order structure functions 
as in Cerutti et al. (2000) and Lindborg (1999). It follows that the Kolmogorov length 
scale (77 — {v^/e)'^/^) was 0.08 mm, the Taylor micro-scale (A) was 2.9 mm, and the 
Reynolds number based on Taylor micro-scale {Re\) was 350. The longitudinal integral 
scale obtained by integrating up to the first zero crossing of the ui correlation function 
was ill ~ 0.091 m. Profiles of mean velocity, rms velocities, Reynolds shear stress, rms 
temperature and heat fiux distributions across the wake agreed quite well with results in 
the literature (e.g., Matsumura & Antonia 1993 and Kiya & Matsumura 1988). 

In the present study, to separate between large and small scales, the box filter is 
applied to the streamwise and cross-wake directions, and the trapezoidal rule is used 
for the spatial integrations. The filtering process consists in a discrete approximation 
to a two-dimensional box filter. In the X2 direction, a four point discretization is used 
for evaluating the SGS fluxes while a three-point approximation is used for the filtered 
derivatives. Filtered velocity and scalar gradients in the X2 direction are evaluated using 
first-order finite differences over a distance h. In the streamwise direction, the box filter 
is approximated using A/s/(wi) sampling points and the xi derivatives are evaluated 
using finite differences over a distance h. The filtering and error analysis is documented 
in Cerutti & Meneveau (2000) and Cerutti et al. (2000). 

Figure 1(a) shows a comparison between the longitudinal spectrum Eii{ki) of the ui 
component and the longitudinal spectrum £^22(^*1) of the U2 component multiplied by 
3/4 at the centerline. Here, ki is the longitudinal wave number. The three vertical dashed 
lines correspond to the filter sizes of A/ry =125 (10 mm), 250 (20 mm), 500 (40 mm). All 
the filter sizes are in the inertial range. Since the noise peak in the longitudinal spectrum 
of Ml is in the far dissipation region, quite removed from any of the filter frequencies 
and scales of interest in this study, no effort is made to remove the noise by additional 
filtering (various attempts such as notch filtering showed no effect on the results) . It can 
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Figure 1. Velocity and temperature spectra at the centerline, in Kolmogorov units: (a) the 
longitudinal spectrum of the iti component (solid line) and the longitudinal spectrum of the 

U2 component multiplied by 3/4 (dashed line), as a function of the longitudinal wavonumbcr, 
Ki. The spectrum of u-z component is below that of ui at low wave rmmbor {nirj < 0.002) 
with a Stroulial peak (Kir/ ~ 0.003), and higher at high wave rmmber (0.1 < Kir; < 0.4); 
(b) the longitudinal spectrum of the temperature (solid line) and of the filtered temperature, 
filtered at A/r? = 125 (dashed line). The inserts in (a) and (b) show the compensated spectra, 

Sli(Kl) 



-2/3, 



l^^, jE22{ki) e~^/^Ki^^ and Eee{Ki) tg^e^^^nX^^ . The three vertical dotted lines 
in (a) and (b) are the wavenumbers corresponding to the filter sizes of A/r? =125 (10 mm), 250 
(20 mm), 500 (40 mm). The straight solid lines are the universal spectra (see text). 



be clearly observed that Eii{ki) « |i?22('*i) as required by isotropy in the inertial range, 
over about one decade of wave numbers. The peak in -E22('*i) kit/ ^ 0.003 is due to 
the periodic von Karman vortex street behind the cylinder. The frequency of the vortices 
is 78.1 Hz, and this gives Strouhal number Si (= fD/Uoo) = 0.211. The Kolmogorov 
constant ck in Eii{ki) = ^ck^^^'^Ki'^^^ is obtained as ck = f|(0.56) = 1.71 in the 
insert. This value is quite close to the standard value of 1.6 (see review by Sreenivasan 
1995) and a similar value of 1.7 was observed by O'Neil & Meneveau (1997). Figure 1(b) 
shows the longitudinal spectrum of the temperature E0g{Ki). As seen in the premultiplied 
spectrum shown in the insert, there is a fairly clear inertial range with a —5/3 slope. This 
slope is steeper than those found in the round jet by Tong & Warhaft (1995) and Miller 
& Dimotakis (1996), and in several other shear flows reviewed in Sreenivasan (1996). It 
is closer to results quoted in Antonia & Pearson (1997), who report a scaling exponent 
of 0.65-0.66 for the 2"''^-order temperature structure function in the heated cylinder wake 
at Rex = 230, or to the spectra for grid turbulence of Mydlarski & Warhaft (1998b). 

The coefficient cq in Eee{Ki) = cgeee~^^^K'^^^^ deduced from this spectrum is about 
1.4. It is significantly higher than results quoted in Sreenivasan (1996) for high Reynolds 
numbers, but is within the range of results from various shear flow measurements reported 
in Antonia et al. 1999. As discussed in Sreenivasan (1996), the dissipation measurements 
are based on assuming scalar isotropy and are thus subject to considerable uncertainty. 
On the other hand, at the centerline we find good isotropy (see section 3) which would 
seem to support our current estimates of the dissipation and cg. While the universality 
of Cg and spectral exponent is not the main subject of this paper, the scatter of re- 
sults certainly supports the view that the scalar spectrum for shear flows at moderate 
Reynolds numbers (e.g. Rex < 1000) depends upon details of the generation of the flow 
(Sreenivasan, 1996). Specifically, it seems that the dependence of the spectral exponent 
and prefactor on Rex is not universal. 

The longitudinal spectrum of the filtered temperature is shown as the dashed curve in 
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Figure 2. Isotropy ratios of the filtered gradient fields and of the mean SGS stresses across 
the heated wake at x\/D — 25 with the filter size of A/rj = 125. The symbols represent: O, 
{522&2)/{SiiSii); □, (Si25i2>/{|5ii5ii); A, {G2G2) / {G^Gi); {T22) /{tu}; M, -{ri2)/(rii). 

figure 1(b). The lobes at scales below the filter size are due to the streamwise box filter 
used in the present study. Spectra of filtered velocity (not shown) have similar shape. 



3. Isotropy of real SGS dissipation and model predictions 

In discussing isotropy of various tensors, we distinguish second-rank and fourth-rank 
tensors. Examples of second-rank tensors are the mean SGS stress (rjj), the filtered 
scalar-gradient (GiGj), and the SGS scalar- variance dissipation —{qiGj). Their isotropic 
form can be written as (ry) = {Tii)Sij, {GiGj) ~ {GiGi)5ij, and (qiGj) = {qiGi)Sij 
respectively. On the other hand, the strain-rate-product tensor {SijSpq) and SGS dissi- 
pation tensor {rijSpq) are fourth-rank tensors. Their isotropic form can be written as 
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SijSpq — — {SipSjq + SiqSjp) 



(1) 



and a similar expression holds for (SijSpg) by replacing r,;j with Sij. To derive the above 
expression, we have used the tensor symmetry (in i — j or p — q) and the divergence-free 
condition {Skk = 0). It follows that {S22S22} = {SiiSn} and (5125*12) = |(5ii5ii). 

Figure 2 shows isotropy ratios of the filtered gradient fields and of the mean SGS 
stresses, for A/rj = 125. In isotropic conditions, all ratios have to be on the unit line 
except (Ti2)/(rii) which should tend to zero with decreasing filter scale since the fraction 
of mean shear stress carried by the SGS scales is expected to vanish when A/£ 0. As 
can be seen, both the SGS stress and strain-rate fields are quite isotropic across the 
heated wake flow at this scale. For the larger filter sizes oi A/rj = 250 and 500 (not 
shown), the trends are similar to figure 2 but slightly less isotropic. 

Next, results are given for the SGS dissipation tensors across the wake. We define the 
following 'isotropy ratios': 



^T22S22j \Ti2S12j \q2G2 ^ 

(rii^ii) f (ni^n) (51G1) 

Figures 3(a) and (b) show the spatial distributions of the SGS kinetic energy dis- 
sipation isotropy ratios /„i2 and the SGS scalar- variance dissipation isotropy ratios Ig, 
respectively, across the wake flow with the different fllter sizes. All moments in this study 
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Figure 3. Isotropy ratios in the SGS dissipations: (a) /«i2; (b) 1$ (see Eq. 2). The symbols 
represent: O, A/r] = 125; □, A/r] = 250; A, A/rj = 500. For clarity, 4*''-order polynomial fits of 
logi vs X2/i are added as solid lines. 



are statistically well converged. For example, running averages such as n ^ X]fc=i li^j 
did not differ from their global average by more than 6% over the last 80% of data (i.e. 
for n G [0.2N,N] where N = 2A ■ 10^). For different filter sizes considered, the SGS 
kinetic energy and scalar- variance dissipations are isotropic only towards the centerline 
of the wake. Towards the edges of the wake, large anisotropics persist (even though there 
the mean shear and mean temperature gradients vanish also). The isotropy ratio of 7„i2 
at a fixed X2/I in fignrc 3(a) clearly increases with the filter size. However, the isotropy 
ratio of the scalar- variance dissipation Ig shows almost no variation with filter size for 
1 < X2/i < 2, where the scalar gradient in the cross- wake direction, is large. 

The increased anisotropy in the outer regions of the wake (see figures 3(a) and (b)) 
is another striking result. The anisotropy may be associated with the highly intermit- 
tent character of the flow there. One may wonder whether the anisotropy arises from 
a superposition of distinct behaviors in the turbulent and the non-turbulent (outer) re- 
gions. Using conditional averaging, O'Neil & Meneveau (1997) already showed that the 
conditional averages of the SGS dissipations in the non-turbulent regions were negligible 
compared to those in the turbulent regions. The implication was that the global aver- 
ages of dissipation could be explained entirely by their conditional mean value inside the 
turbulent part: {TijSij) ~ T{TijSij)T, where T{x2) is the intermittency function (i.e. the 
fraction of time the signal is turbulent at any given X2), and the subscript "T" stands for 
avc;raging conditioned on 'turbulence' (sec; O'Neil & Meneveau 1997 for details). When 
replacing the conditional averages multiplied by F in the expression for isotropy ra- 
tios, r cancels from both numerator and denominator. This behavior suggests that the 
anisotropy ratios shown in figure 3 are equal to those inside the turbulent regions alone, 
and that their rise in the outer parts of the wake cannot be explained by contributions 
coming from the non- turbulent parts. However, it is still possible that non-trivial con- 
tributions to the scalar-variance dissipations could originate at the interface separating 
turbulent from non-turbulent regions. 

To highlight the variations with filter scale more directly, in figure 4 we plot the 
isotropy ratios as a function of scale, for the transverse location X2/£ — 1.44 close to the 
peak mean temperature gradient. As is evident, the isotropy ratios of the kinetic energy 
dissipation decrease towards unity as the filter size decreases, whereas the isotropy ratio 
of the scalar-variance dissipation remains unchanged near Ig ^ 2 as the filter size is 
decreased. As seen in figure 3, near the centerline, where the mean shear and scalar 
gradient vanish, the SGS isotropy ratios are all near unity. 
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Figure 4. Isotropy ratios of kinetic energy and scalar- variance SGS dissipations at x-z/l = 1.44, 
where the mean temperature gradient peaks. The open symbols (and dashed lines) represent 
isotropy ratios for the SGS dissipation of kinetic energy, namely O, 7^22; and □, /ui2; the closed 
triangles ▲ (and solid line) are for the scalar dissipation I0. 



Therefore, wc conclude that in terms of the most important of the resolved-SGS scale 
interactions (the mean SGS dissipation), the scalar field maintains strong anisotropy 
in the presence of a mean scalar gradient. Conversely, the velocity field has the trends 
expected from approach to isotropy at small scales. 

Next, we quantify the isotropy level of model predictions, when the data are analyzed 
in an 'a-priori' sense, i.e. by replacing Tjj and Qj above by model expressions. First, the 
standard eddy-diffusion model is considered, i.e.. 



Smag 
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Smag 



2Pr-'{Cs/^f\S\G^, (3) 



where |^| = {2SmnSmn) ' is the modulus of the resolved strain rate, Cs is the Smagorin- 
sky coefficient, and Prggs is the SGS Prandtl number. Consequently, and independent of 
the model coefficients, the isotropy ratios from the eddy-diffusion model are: 



rSmag 
^u22 



{\s\s.. 



22'->22 



\S\SnSn 



rSmag 
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(4) 



In computing the filtered strain-rate magnitude from the data, the following approxima- 
tion is used (this is a 2D extension of the ID approach of O'Neil & Meneveau 1997, and 
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was also used in Liu, Katz & Meneveau 1999): l^l 

This approximation is not expected to affect the accuracy of our measured isotropy ratios 
significantly. The reason is that l^j and its 2D approximation are scalars which multiply 
equally all terms of the squared velocity and scalar gradients and is further supported 
by the isotropic behavior of second-order moments of Sjj shown in figure 2. The isotropy 
ratios of the SGS dissipations from the eddy-diffusion model are shown in figure 5(a). 
Results are near unity almost independently of scale, including the passive scalar dissi- 
pation. Profiles across the wake (not shown) also arc near unity. This result is consistent 
with the observed isotropy of the square gradient tensors shown in figure 2 (the only 
difference here is the additional |5| factor). Hence, using the eddy-diffusion model one 
would (incorrectly) predict SGS isotropy since the resolved second-order moments are 
isotropic. 

The second modeling class to be considered here is the 'nonlinear model' or 'tensor 
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Figure 5. Isotropy ratios of components in SGS dissipations from models at X2/£ = 1.44, as a 



function of filter scale: (a) eddy-diffusion models: O, 7^ 
nonlinear models: O 
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eddy-diffusion model' (Leonard 1974, Clark, Ferziger & Reynolds 1979, Liu, Meneveau 
& Katz 1994, Leonard 1997, Borue & Orszag 1998 and Meneveau & Katz 2000). The 
nonlinear model reads as follows: 



if 



Cnl 



^ 9 do du-j 



(5) 



dxk dxk ' " * dxk dxk ' 

where C„; and Cni^ are the nonlinear model coefficients. Therefore, the corresponding 
isotropy ratios from the nonlinear model are defined as follows: 



rnl = 
-'u22 — 



dU2 dU2 Q 

dui dui Q 



/ dui 

-'«12 = TT^ 



du2 q 
k dx ^ 



12 



dui dui c 
dxk dxk "^ll 



7-nl _ \ ^^fe ^^fc 
ia = 



06 dui (J 

dxk dxk ^ 



(6) 



where k varies from 1 to 2 in the analysis of our 2D data. 

The isotropy ratios of the SGS dissipations from the nonlinear model are shown in 
figure 5(b), as a function of filter scale. As is apparent in comparing with figure 4, the 
main features of anisotropy and filter-size dependence are correctly reproduced qualita- 
tively. Quantitatively, the levels of anisotropy are overestimated. For instance, the level of 
anisotropy for the modeled SGS scalar- variance dissipation appears to stay near Ig^ ~ 2.5 
as opposed to ^ 2 for the real SGS scalar dissipation. 



4. Summary and conclusions 

In studying passive scalar statistics in a turbulent shear flow with a mean temperature 
gradient we focus on statistics of interest to subgrid modeling and large-eddy simulation. 
In order to obtain the filtered and subgrid velocities and temperatures, a probe array 
composed of four X-wire and four cold-wire sensors is used and two-dimensional box 
filtering in the streamwise and cross- wake directions is applied to the data. The isotropy 
ratios of the SGS kinetic energy and scalar- variance dissipations are investigated as func- 
tion of position in the flow and as a function of fllter scale. Both dissipations are isotropic 
independent of filter size near the centerline where there is no mean shear or scalar gra- 
dient. However, at locations with high gradient (and also further towards the outer wake 
regions), we find that the scalar- variance dissipation remains highly anisotropic, indepen- 
dent of filter size. Conversely, the kinetic energy dissipation tensor approaches isotropy 
as the filter size is decreased. A mechanistic explanation of the observed trends in terms 
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of possible orientations of ramp and cliff structures is not evident to us at this time. 
The persistence of scalar anisotropy even at small scales is consistent with prior results 
for structure functions and gradient statistics of unfiltered turbulence (Mestayer 1982, 
Sreenivasan 1991, Warhaft 2000). Present results quantify the impact of this anisotropy 
on the interactions among large and small scales in the context of SGS modeling and 
LES. 

We find that the predictions of eddy-diffusion models arc much more isotropic than the 
real phenomenon. This may at first glance seem obvious since eddy-diffusion models are 
often invoked under the banner of isotropy. However, the present result does not arise from 
an explicit isotropy assumption built into the model, but because the resolved gradients 
have second-order statistics that are isotropic. For instance, if the scalar gradient tensor 
(GiGj) had been found to be anisotropic, it would have implied anisotropic behavior of 
the eddy-diffusion model's predictions. Instead, the isotropy that exists in the filtered 
scalar gradients is incorrectly applied to model the anisotropic statistics of the SGS heat 
flux. The main problem for the eddy-diffusion model appears to be that it uses second- 
order statistics to model third-order statistics. Conversely, the anisotropy exists in the 
third-order moments that arise from the velocity and scalar product (gj) multiplied by the 
scalar gradient (Gj), but is not discernible in the second-order statistics of Gj alone (even 
when modulated by the strain-rate magnitude). The anisotropy is clearly discernible, 
however, in the third-order moments consisting of the filtered velocity gradients times 
scalar gradients squared that arise in the expression for modeled SGS dissipation of scalar- 
variance using the nonlinear model. These expressions are able to reproduce the detailed 
phase relationships among the velocity and scalar field that govern the SGS dissipation 
(cascade) of scalar variance. We remark that the ovcrprediction of anisotropy by the 
nonlinear model (and its underprediction by the eddy-diff\ision model) is reminiscent 
of the opposing trends of these two models documented in Liu et al. (1999) for rapidly 
strained turbulence in cold-flow. The opposing trends suggest that a linear combination 
of the two models, i.e. the 'mixed model', can be tuned to reproduce the correct amount 
of anisotropy (for a discussion of the application of mixed models in LES, see Meneveau 
& Katz 2000). However, note that the data for /„i2 and /„22 at A/r; = 250 and 125 do 
not provide clear justification for choosing the mixed over the non-linear model. 

Finally, it is stressed that current results are obtained in a single fiow for a single 
moderate Reynolds number. Even if the present evidence in figure 4 of an essentially 
scale-independent anisotropy for the scalar dissipation appears to be quite strong, the 
results could change in another fiow, or at higher Reynolds numbers (we recall that 
according to Sreenivasan 1996, universal behavior for the scalar requires Re\ above 1000 
or so). These considerations serve as motivation for further work in this area. 

We thank Professors Z. Warhaft and L. Mydlarski for useful comments and suggestions 
about the cold-wire calibration. This work was financially supported by the National 
Science Foundation (grant CTS-9803385). 
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